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Abstract 

We investigate the core structure of radiatively cooling intracluster gas, using a hydrodynamics code. 
We calculate evolution of model clusters of the initial core radii 160-300 kpc until the initial central cooling 
time, and analyze the resultant clusters using the double /3-model as done by observational studies. It is 
found that the core-size distribution thus obtained shows two peaks ~ 30-100 kpc and ~ 100-300 kpc and 
marginally can reproduce the observed distribution which exhibits two distinct peaks around '^ 50 kpc and 
~ 200 kpc. This result may suggest radiative-cooling origin for small cores, while cooling is yet insignificant 
in the clusters of large cores. It should be noted that the small core peak is reproduced by clusters that 
are still keeping quasi-hydrostatic balance before the initial central cooling time has elapsed. 
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1. Introduction 

Profiles of X-ray emitting hot gas in clusters of galaxies 
have been studied often by using the so-called /3-model, 
which is an isothermal hydrostatic gas model consisting 
of a core and envelope. From the /3-model analyses of 121 
clusters including nearby clusters (Mohr et al. 1999), Ota 
and Mitsuda (2002) obtained an interesting result that 
the distribution of the core radii of the intracluster gas 
exhibits two distinct peaks at '^ 50 kpc and ^ 200 kpc 
for Ho = 70 km s'^ Mpc"^ (see also Ota, Mitsuda 2004). 
Akahori and Masai (2005) (hereafter AM05) investigated 
correlations of the core radii with various properties of the 
clusters, and found that the radii in the larger core group 
around the peak '^ 200 kpc are marginally proportional 
to the virial radii of the clusters and therefore the origin 
may be attributed to simple self-similar collapse. On the 
other hand, in the smaller core group around ~ 50 kpc, no 
clear correlation is found between the core and virial radii, 
suggesting some other origin of the small core formation. 

AM05 investigated several possibilities of the origin 
of the small cores, not only by examining correlations 
among the observational quantities but by simulating the 
/3-model with a hydrodynamics code. Many of the small 
core clusters possess central cD's or giant ellipticals, but 
the simulation shows that the gas core under their pres- 
ence is ~ 40 kpc at most, and is too small to account for 
the observed range ^ 40-80 kpc. AM05 suggested another 
possibility: the effects of radiative cooling on the core size. 
Radiative cooling time is shorter in the core than in the 
ambient region, because the density is higher. As the core 
cools, the ambient gas then could inflow to compensate the 
pressure loss inside. Although the ASCA. Chandra and 
XMM-Newton observations (e.g., Makishima et al. 2001; 
Lewis et al. 2002; Peterson et al. 2001) suggested much 
smaller amount of cooled mass than expected for the clas- 



sical cooling flow model (see Fabian 1994 for a review), 
this process increases the gas density toward the cluster 
center and is likely responsible for the small cores ob- 
served. AM05 found that the central cooling time, icooi; 
is significantly shorter than the Hubble time for 48 out of 
50 small core clusters. 

Thermal evolution of clusters has been studied by many 
authors. Most of the papers addressed processes that 
would slow or inhibit the onset of radiative cooling in- 
stability, such as electron conduction and activities of ac- 
tive galactic nuclei. Until an elapsed time comparcble to 
tcooi, however, the gas appears to be cooling rather slowly 
with the temperature ~ keV at the cluster center. For 
example, Ruszkowski and Begelman (2002) studied the 
heating, conduction, and minimum temperature in cool- 
ing flows with a simple spherical model and showed that 
line and free-free cooling in the cluster center leads to slow 
cooling until the initial central cooling time has elapsed. 
At this stage feedback is not yet very important. Masai 
and Kitayama (2004) showed that when the intracluster 
gas is cooling with quasi-hydrostatic balancing, the tem- 
perature of the cooling gas appears to approach a constant 
value toward the cluster center. 

In the present paper, we investigate thermal evolution 
of clusters, paying attention to the core radius of the intr- 
acluster gas undergoing radiative cooling. Our interest is 
the origin of the small cores or the observed core-size dis- 
tribution which exhibits two distinct peaks. We simulate 
cooling clusters using a hydrodynamics code (AM05), and 
analyze the core radii applying the /3-model, on which the 
observed results are based. In section 2 we describe the 
model. In section 3, we present the results of calculations 
and discuss the properties of the cooling cores. We give 
some concluding remarks in section 4. 



Model and Calculations 
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Table 1. Best-fit parameters for the King profile. 



We calculate the evolution of radiatively cooling intra- 
cluster gas "which is initially in hydrostatic balance under 
the presence of dark matter. We analyze the resultant 
clusters using the /3-niodel in the same manner as the ob- 
servations (Ota, Mitsuda 2002): the core-radius distribu- 
tion having two distinct peaks was obtained by analyses 
using the /3-model. Since the /?-modcl is based originally 
on the King model of the coUisionless matter (see e.g., 
Sarazin 1986), we apply primarily the King model for the 
dark matter and galaxies of a cluster. We also mention 
the results in the case of the NFW dark matter model 
(Navarro et al. 1996), because such a cuspy profile has 
been suggested for dark halo mergers or cluster formation 
by numerical simulations. As for the thermal evolution of 
the intracluster gas, the two dark matter models show the 
similar results to each other. In the following we describe 
our models and calculations along the King dark matter 
model. 

With the gravitational potential $ of the total mass 
within radius r, gravitational balance of a spherically sym- 
metric cluster is written as, 
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where "dust" means coUisionless components, i.e. dark 
matter and galaxies, k is the Boltzmann constant, m the 
proton mass, fi the mean molecular weight for which we 
take fj, = 0.6, p the mass density for each component. T 
the gas temperature, and ct^ is the radial velocity dis- 
persion which is equal to the line-of-sight velocity disper- 
sion for the isotropic dust. If isothermal distribution is 
attained for every component, i.e. dlnT/dhir ~ and 
dlna'l^/dlnr ~ 0, 
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follows equation (1). 

For comparison with the observations of core sizes, in 
modelling clusters wc apply the density profile function 
based on the /3-model (Cavaliere, Fusco-Femiano 1976), 



p(r)=po[l + (r/rc) 



21-3/3/2 



(3) 



where po is the central density, Vc is the core radius, and 
/3 is a parameter to represent the envelope slope. For the 
consistency with the /3-model for the gas, we adopt the 
King model for the coUisionless components. Although 
the original King model (King 1966) is not of simple ana- 
lytic form, the profile may be approximated with equa- 
tion (3), which is sometimes called approximate King 
model (see e.g., Sarazin 1986). Table 1 shows the best-fit 
parameters for the King profile with equation (3) in the 
range r/rx = 0.0-10.0, where tk is the core radius of the 
King model. 
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Fig. 1. Radial density profiles. The solid and dashed lines 

2/3 
represent p» oc pK and pg oc Pj^ , respectively. The dotted line 

is the approximate King profile (equation (3) with /3 = 1). 



We investigate four models of clusters which consist of 
dark matter, galaxies, and gas with their mass ratio Mgai : 
Mgas : Mdm = 1 : 5 : 30. The galaxies and dark matter are 
represented by their gravitational potential components 
in the momentum equation (equation (16) in AM05). We 
employ the original King model for their density distribu- 
tions, and then pg oc p^.^"^"" with /3spoc(— /3prof) = 2/3 for 
the initial gas profile (equation (2)). The initial profiles 
are shown in figure 1, where the tidal radius rt = rvir, i.e. 
P:,{r > Tvir) = 0, and the core radius rx — J'c- The core 
radius of the gas well reflects the gravitational potential 
consistently in the context of the /?- model (AM05). 

The four model-clusters have different initial core radii 
of the gas, re = 160, 200, 250, and 300 kpc, and sat- 
isfy the self-similar relation with their virial radii, as 
x-vir = rvir/rc = 15, or r^i, = 2.40, 3.00, 3.75, and 4.50 Mpc 
(listed in table 2). These values are on the dashed line in 
figure 5 of AM05, which represents the r^i^/rc relation ob- 
tained for the large core group peaked at ~ 200 kpc in the 
observed core-size distribution. The initial gas tempera- 
ture (isothermal), T, under the hydrostatic equilibrium is 
higher than the virial temperature, Tvir, as (3spccT ~ Tvir 
for xl-Jil -K xl-^) ~ 0.996 (see AM05). 

We carry out hydrodynamical simulations using 
a Smoothed Particle Hydrodynamics (SPH) code 
(Hernquist, Katz 1989; Monaghan 1992) with 95,000 bod- 
ies. The smoothing length is estimated to be ft. = 17.1 kpc 
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Core Structure of Intracluster Gas 

Table 2. Initial properties of the simulated clusters. 



/J-model parameters** 

rvir Mvir* T re Ugf) (3 icool** 

Model (Mpc) (lOi^Afo) (keV) (kpc) (O.Olcm^^) (Qy^.) 

(a) 2.40 0.555 3.12 171 2.15 0.83 0.689 

(b) 3.00 1.08 4.87 209 2.16 0.82 0.800 

(c) 3.75 2.62 7.61 260 2.10 0.81 1.086 

(d) 4.50 3.66 11.0 310 2.04 0.80 1.303 

* The virial mass. 

** Values are estimated after 1.6 Gyr of the hydrostatic balance test without cooling. 



and 9.97 kpc for the core region (rig ^ 0.02 cni^'^) and 
cooled cluster center (rig ^ 0.1 cm"'^). respectively, and 
sufficiently small compared with the core size concerned 
here. Here, Tig = pg/^m is the number density of the 
gas. For the general artificial viscosity we adopt the 
Monaghan-Gingold value (Monaghan. Gingold 1983; see 
also Hernquist, Katz 1989) with the coefficients (a, (3) 
= (1.0, 2.0), and for the gravitational softening we take 
V$ oc l/(r2 + e^) with e = O.l/i. 

Radiative cooling is taken into account for the energy 
equation as 

dui Pi 



dt 



rV • u, - PiA{T,] 



(4) 



where Ui is the internal energy per unit mass, Vi is the 
velocity, and A(Ti) is the cooling function for the «-th 
SPH-particle. No additional heating processes are con- 
sidered because the evolution of cooling cores is of our 
interest. Neither is electron conduction, which is not ef- 
fective within cooling cores. We use the cooling func- 
tion of Sutherland and Dopita (1993), which includes line 
emission with the metalicity Z~0.3Zq, by approximation 
as A{Ti) = 3.0 X 10"27y^i/2 ^^^ ^^s g-i ^^ ^^le tempera- 
ture range ~1.5-10 keV. Actually, until t ~ tcooi of prac- 
tical interest, the temperatures of the cooling gas are still 
above ^ 1.5 keV where bremsstrahlung dominates. Here, 
icooi = 3ngoA:T/(ngo^A) is the initial cooling time at the 
cluster center. 

Before going to calculations of evolution with radia- 
tive cooling, we examine the hydrostatic balance of the 
gas with A{Ti) = in equation (4), as shown in figure 2. 
Gravity on the gas and the pressure gradient balance with 
each other, and the initial gas profile is kept at least until 
3.2 Gyr, which is enough long compared with the dynami- 
cal (free-fall) timescale, t^ = (32/37rGp)^/^ '^1.6 Gyr at the 
cluster center {ugo = 0.02 cm^'^). The gas is kept nearly 
isothermal so that the profile is well represented by the 
/5-model, although the temperature declines somewhat at 
r > bvc due likely to adiabatic expansion of the outermost 
envelope. At r < 5rc of practical interest, the hydrostatic 
gas is represented well by the /3-model as follows. From 
equations (2) and (3), the /3-model should give a relation 
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Fig. 2. Initial hydrostatic balance of the simulated clusters. 
We obtain the clusters after 1.6 Gyr of the hydrostatic balance 
test without cooling. For each box, the upper panel represents 
the gravity on the gas and the pressure gradient (the thin 
and thick line, respectively, in a unit 1.79 X 10^® g cm s~^), 
and the lower panel represents the density (the thin line, 
0.01 cm~^) and temperature (the thick line, keV). 



or /3 ~ /5*/3spcc. Actually, /3 - 0.80-0.83 (table 2) is 
consistent with /?*/3spcc ^ 0.8 for /3* ~ 1.2 (table 1) and 

Pspcc — ^/ <^- 

3. Results and Discussion 

3.1. Evolutions of the (3 -Model Parameters 

We apply the double /3-model as well as the single (3- 
model to simulated clusters. The double /3-model is the 
superposition of two single /3-nLodels as 

pir) = pi[l + (r/ri)2]-3/3i/2 ^ „.jt ^ rr/r-.^2i -3/32/2 



-P2[l + {r/r2y 



(6) 



/3prof = P/(3. ^ A 



spec 5 



(5) 



While in observations the parameters are inferred from the 
surface brightness profile, which depends substantially on 
the density profile as oc p^T^/^, in our calculations they 
are obtained straightforwardly from the density profile. 
We obtain the best-fit parameters of the "outer" com- 
ponents (pi, ri, and /?i) using the data in the range 
1.0 < r/rK < 5.0, and then obtain the "inner" ones (p2 
and r2), assuming P2 = Pi as done by Ota and Mitsuda 
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Fig. 3. Thermal evolutions of the /3-irLodcl parameters: slope 
parameter (the top panel), number density of center (the mid- 
dle, in a unit 0.01 cni~^), and core radius (the bottom, Mpc). 
The thin and thick solid lines represent the outer and inner 
components of the double /3-model, respectively. The dotted 
lines represent the components of the single /3-model. The 
arrows represent the cooling time of each cluster. 

(2002; 2004), using the data of 0.0 < r/rx < 5.0 including 
the outer component with pi, ri and /3i fixed. 

Figure 3 shows the time evolution of the parameters for 
the single and double /3-models after the gas starts to cool. 
We see in the bottom panel that the outer core radii (the 
thin solid lines) decrease very slowly compared with the 
core radii of the single /3-model (the dotted lines). This 
implies that the outer component keeps the initial value 
since thermal evolution of the gas is predominant in the 
inner component. Actually, /9i, which is determined by 
the outer component, is roughly kept the initial value, 
although (3 of the single /?-model decreases by about 10% 
within fcooi (the top panel). The behavior is also seen 
in the central density. While the increase in the central 
density of the outer component is as small as a factor of 
~ 2, the central density with the single /3-model increases 
exponentially (the middle panel). 

From the time evolution of the core radii, we calculate 
the time during which a cluster would have the core ra- 
dius between r and r -\- Ar and estimate the population or 
relative number of clusters that would fall into a certain 
range of core radius. The core size distribution thus ob- 
tained for the outer and inner cores is shown in figure 4. 
Our calculations reproduce the observed core-size distri- 
bution or two distinct peaks except for the details such 




Fig. 4. Distribution of core radii. The thick solid lines rep- 
resent the relative probabilities within tcool of the four sim- 
ulated clusters, estimated for the outer and inner cores with 
the bin of Ar = 20 kpc and 10 kpc, respectively. The dotted 
line represents the distribution of core radii of 121 clusters 
(same as figure 4 of AMOS). The thin solid lines represent the 
relative probabilities for the case of the NFW dark matter 
potential. 

as their widths or the tails. This implies that radiative 
cooling is a possible cause of small cores of clusters, while 
cooling is yet minor in large core clusters. It is interesting 
that even the small core peak is produced by the clus- 
ters of t ^ icooi, i-C. moderately cooled clusters that are 
keeping quasi-hydrostatic balance. 

In the observed distribution of figure 4 one may notice 
four clusters which have very small cores. It is unlikely 
that such small cores as r^ ^ 30 kpc are reproduced by 
radiative cooling within t ^ tcooi- Our cooling function 
underestimates by ~ 40 % at T = 1.5 keV compared with 
the function of Sutherland and Dopita (1993). If applied 
their original function, the small core peak might have a 
tail toward the small core end. However, this is unlikely 
because such a very small core is quite transient (figure 3). 
This may be related with the fact that clusters of i ^ tcooi 
are out of quasi-hydrostatic balance. Another possible ex- 
planation of the very small cores is the presence of central 
gravitational source; the four clusters evidently have cen- 
tral cD's or giant ellipticals. The presence of such galaxies 
is likely responsible for the very small cores (sec AM05). 

The lack of clusters of core radii ~ 120 kpc is seen in 
the calculated distribution as well. Although the outer 
cores of ^ 100-300 kpc of simulated clusters can explain 
the population around the peak of the observed large core 
clusters, the latter exhibits larger core tail extended to 
Tc ~ 1 Mpc. If we simulate clusters which have initial 
huge cores ^ 0.4 Mpc, we might get the cooled cores 
of ~ 120 kpc following the self-similarity (see below). 
However, if such huge core clusters are in hydrostatic equi- 
libria as represented by the /?-inodel, the clusters would 
have unlikely virial radii as large as r^-„ ^ 6 Mpc. In fact, 
most of the huge core clusters deviate from the self-similar 
relation, as shown in figure 5 of AM05. 



No. 



Core Structure of Intracluster Gas 



A considerable factor for the origin of such huge cores 
is mergers. When merger takes place, the gas could form 
a large flat core with a steep envelope (large /3) . Actually, 
in AMOS sample, 6 out of 7 clusters of /3 > 1 have cores 
of Tc > 0.6 Mpc, and the average of /3 in the 17 clusters 
(including 14 irregulars) of r^ > 0.4 Mpc is 0.98, which is 
significantly greater than the average 0.65 in 121 clusters. 

We find in figure 3 that the single /3-model represents 
the intermediate profile between the inner and outer pro- 
files of the double /3-model. If cooling is minor and the 
self-similarity remains, the profile would be represented 
by the single /3-model of a large core. With cooling the 
core, but yet minor in the outer, the profile is being that 
represented better but transiently by the double /3-model, 
and eventually comes to that represented by the single 
/3-model of a small core. As a result, the core-size dis- 
tribution exhibits two distinct peaks along the thermal 
evolution. Two distinct peaks are not seen for the cluster 
sample of rcdshift between 0.4 and 1.3 (Ettori et al. 2004). 
The reason may be that their analysis was done with the 
single /3-model; otherwise, the cooling might be yet minor 
in such very distant clusters. We confirm that two distinct 
peaks do not appear if applied the single /3-model to the 
simulated clusters. The estimation of the inner core ra- 
dius depends on /3i obtained from the outer slope. If /3i is 
taken to be 10 % smaller than that in the present calcula- 
tions, the resultant clusters might have 10 % smaller cool- 
ing cores, as expected from the /3-rc correlation described 
in figure 1. Thus, the small-core distribution merely shifts 
toward the smaller side. 

We also investigate clusters of the NFW dust profile 
(Navarro ct al. 1996), where the parameters of the NFW 
profile arc determined so that the clusters have nearly the 
same initial density and temperature profiles of the gas, 
and analyze the core sizes in the same manner as above. 
The obtained core-radius distribution is similar except for 
a few % tail extended to ^ 20 kpc (figure 4). There is no 
appreciable difference between the NFW and King dark 
matter models in the cooling gas distribution. 

The evolutions of the density and temperature profiles 
(figure 5) imply the self-similar evolution in the core re- 
gion, since the evolution is predominant in the inner com- 
ponent. We try to approximate the evolution of inner- 
core radii with a parameter, e, and the normalized time, 

T = l/Ccool? <iS 



1^2 = ''2,0(1 — <^t). 



Figure 6 shows that equation (7) can represent the inner- 
core evolution, i.e. thermal evolution of the small core 
can be described by the self-similar relation, normalized 
by icooi and r2,o which are obtained from the initial outer 
component. This self-similarity may be seen also in the 
clusters having different initial central densities, because 
in the /3-model the central density determines the ampli- 
tude of the profile while the flat-core region clearly exists; 
the core radius is related mainly to the envelope slope (see 
figure 1). 

Such a self-similarity is, however, likely lost in the ob- 
served small cores. At the initial state, the inner-core size 
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Fig. 6. Thermal evolutions of inner-core radii. The dashed 
lines represent equation (7) with e = 0.3, 0.5, 0.7. 

can be represented by the outer-core size, as 

r2,Q^ari,o, (8) 

where a '^ 0.3 is obtained from figure 3. While the thermal 
evolution is predominant in the inner component, outer 
core is roughly constant, i.e. ri ~ ri^Oi so that 6r2 = r2,o — 
r2 is written with equations (7) and (8) as 



6r2 ' 



t-cool 



-t. 



(9) 



From equations (10) and (11) in AMOS, the initial central 
cooling time is given by 

-1/2 



tr. 



CX [13 



1 



ri,o. 



(10) 



Therefore, equation (9) is written as 
Sr2 (X ae(3-^/^t 



(11) 

for the clusters with Xvir ^ 1 (AM05). Equation (11) 
suggests that the self-similar relation in the small cores 
may be lost for not only the different age of the cluster 
from the last major merger, i.e. the time when the cluster 
formed to the current size and down to the hydrostatic 
equilibrium, which is labeled by i, but the various values 
of /3 from ^ 0.6 to ^ 1.0 in observed clusters as mentioned 
in AMOS. 



(■7) 3.2. Evolution of Hydrostatic Structure 



Calculations show that the cooling gas approximately 
keeps the hydrostatic balance between the gravitational 
force and pressure gradient at each step until ~ icooi (fig- 
ure 7). We compare our result with quasi-hydrostatic 
cooling discussed by Masai and Kitayama (2004). They 
suggest that M^ ~ -^r/^cooi at every r or 



C{r) = {Mr/icoo\)/Mr ^ constant, 



(12) 



where Air = 4:Tr pgr'^dr/dt is the continuity equation, /\/,. = 
4:Trr^Pg/3 and icooi = SrigkT / {rig^ A) are the mass of a uni- 
form gas sphere and the cooling time for the local values 



T. Akahori and K. Masai 



[Vol. 




E 1 1 ' ' E 


t^i:^:^^^^,^^^^ ^ 













0.02 0.05 0.1 0.2 
R [Mpc] 



0.02 0.05 0.1 0.2 
R [Mpc] 



0.02 0.05 0.1 0.2 
R [Mpc] 



0.02 0.05 0.1 0.2 
R [Mpc] 



Fig. 5. Thermal evolutions of the density and temperature profiles of the four simulated clusters. The lines are drew for intervals 
of about ~ 200 Myr until tcool • Dotted lines represent the profiles for continuing calculations over tcool > which are excluded for the 
analysis in figure 4. 
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Fig. 7. Hydrostatic balance of cluster (b) after 0.4 (0.5 icool)i 
0.6 (0.75 tcool), 0.8 (tcool), and 1.0 (1.25 tcool) Gyr. The lines 
are the same descriptions as figure 2. The dotted lines are the 
initial values of the density and temperature. 

at radius r, respectively. As shown in figure 8, our calcula- 
tion shows C{r) is nearly constant. Until ~ icooi, constant 
C is seen commonly for all the model clusters, not only 
the King case but also the NFW case. 

Masai and Kitayama (2004) show some properties 
of quasi-hydrostatic cooling by approximating the den- 
sity/temperature profiles in the power-law form. We ex- 
amine our results in the same manner with pg oc r" , T (xr^ 



and inflow velocity Wj, 



r-C 



Since C{r) 



„l+a 



quasi-hydrostatic cooling means 1 + a — ry/2 — ^ • 



-,,/2-C 

or 



-l + ?;/2 + C- 



(13) 
40 kpc: a = -0.42 



Such a relation is actually found at 

is comparable to -1-f ?7/2-f C 0.41 with (77, C)==(0.28, 

0.45). At ~ 70 kpc, however, a = —0.63 is smaller than 
-1 + 77/2 + C- -0.54 with (?7, C)=(0.25, 0.34). Therefore, 
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quasi-hydrostatic balance may be satisfied marginally. 
When the cluster center cools rapidly at t > <cooi (figure 5), 
the balance breaks because a large amount of inflow is re- 
quired to maintain the balance. This also leads to the fact 
that the spatial resolution in calculations becomes prac- 
tically worse than estimated simply from the number of 
SPH particles. 

At r < 100 kpc, the gas inflow velocity Win ^ O.lcs, where 
Cs is the sound speed. The mass inflow rate Mr(100 kpc) ~ 
200 Mq is decreasing to Mr(20 kpc) ~ 10 Mq, which are 
much smaller than expected for the classical cooling flow 
model and consistent with the quasi-hydrostatic cooling 
model by Masai and Kitayama (2004). It should be noted 
that they consider C < for smooth inflow, but we find 
in figure 8 that the inflow velocity gradually decreases 
toward the cluster center, i.e. C > rather than (^ < 0. It is, 
however, consistent with their prediction that icooi oc r^~'» 
is well satisfied in r < 100 kpc. 

C > implies that the gas is heated by viscosity. 
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We confirmed the convergence of the solution by using 
SPH-particles greater than ~ 30,000, which is also suffi- 
cient resolution to remove improperly accelerating cooling 
(Springcl; Hcrnquist 2002). The gas inflow might be ex- 
aggerated by underestimating the viscous heating at the 
dead center of the cluster, though the influence reaches 
up likely to ~ 2h from the center because SPH-particles 
within 2h at a given point contribute to the smoothed es- 
timate in the present work. Therefore, at each step until 
~ icooi, quasi- hydrostatic balance is likely attained with- 
out some heat sources such as electron conduction and/or 
AGN activities. 

4. Concluding Remarks 

We investigate the thermal evolution of cluster cores: 
the core size of the intracluster gas can vary as the gas 
cools radiatively and flows toward the cluster center. In 
order to compare with the observed core radii, we apply 
the /3-model or double /3-model, as done by the obser- 
vational studies, to analyze the gas profile, although the 
model is not always good for cooling cores of more or 
less center-peaked profiles. Thermal evolution of the gas 
may be classified into the following three stages: (I) at 
T = t/tcooi ^ 0.5, the profile is represented by the single 
/3-model with a large core or by the outer-core dominated 
double /3-model, (II) at 0.5 ^ r < I, the profile is repre- 
sented well by the double /?- model, and (III) at t ^ 1 , the 
profile is represented by the single-/3 model with a small 
core or by the inner-core dominated double /3-modcl. 

Until ~ icooi, which is the central cooling time for the 
initial gas profile, the gas cools keeping the hydrostatic 
balance between the gravity on the gas and the gradient of 
the thermal pressure. This evolution is explained by the 
quasi-hydrostatic cooling model proposed by Masai and 
Kitayama (2004) . The properties of the cooling gas in our 
calculations, such as the mass inflow rate, radius depen- 
dence of the local cooling time and constant C, are consis- 
tent with this model in the regime where quasi-hydrostatic 
structure is attained. We find that the inflow velocity of 
the gas decreases toward the cluster center, and its mass 
flow rate is about 10 Mq at 20 kpc. Quasi-hydrostatic 
condition is marginally satisfied up to i ~ icooi , and then 
the cluster center cools rapidly. Regarding the relation 
with the core-size distribution, it may be an important 
clue for understanding thermal properties of clusters that 
even the small core peak is produced by clusters of i ^ tcooi , 
in other words, by clusters that are still keeping quasi- 
hydrostatic balance. 

Analyzing the simulated clusters with the single and 
double /3-modcls, we demonstrate the core size distribu- 
tion exhibits two distinct peaks at ~ 100-300 kpc and at 
~ 30-100 kpc with a valley at ~ 120 kpc until tcooi and is in 
fairly good agreement with the observed distribution (Ota 
and Mitsuda 2002; Akahori and Masai 2005). This imphes 
that the origin of small cores can be explained by radia- 
tive cooling or thermal evolution of the gas. However, a 
question remains: the self-similarity is kept in simulated 
clusters through the thermal evolution with cooling, but 



seems to be lost in the observed small core clusters for 
the various values of (3 from ^ 0.6 to '-^ 1.0 (Akahori and 
Masai 2005). 

The authors would like to thank Tetsu Kitayama for his 
help for computation resources, and him and Naomi Ota 
for their useful discussions. 
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